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ABSTRACT 

If gamma-ray bursts (GRBs) are cosmological they need not necessarily be spa- 
tially coincident with a host galaxy — in many cosmological models the progenitors are 
ejected to large distances from their parent galaxies. While the optical transient of 
GRB 970228 may be close to a faint galaxy, there are several other galaxies in the 
field which could be the host if the progenitors are ejected to ^ bOh^^ kpc. If GRBs 
come from such a widely distributed progenitor population then it should be possible 
to search for a statistical correlation between optical transients and nearby galaxies. 
We provide a statistical framework for quantifying a correlation between potential 
host galaxies and GRBs, taking into account that galaxies are clustered. From simu- 
lations, we estimate that if GRBs occur within I < 100 kpc of their host, then optical 
observations {B < 24.0) of the galaxy field near just 6 optical transient positions are 
required to find a correlation with potential hosts at the 95 percent confidence level. 
The methodology constructed herein is extendible to any cross-correlation analysis 
and is especially useful if either set of objects is clustered. Ultimately, if the distance 
scale to GRBs is established as cosmological, this analysis can be used to determine 
the spatial distribution of GRBs near their host galaxies. 

Key words: gamma-ray bursts — faint galaxies — correlation statistics. 



1 INTRODUCTION 

After nearly thirty years since the discovery of gamma-ray 
bursts (GRBs), a fading afterglow has been detected in X- 
ray and optical wave-bands (Van Paradijs et al. 1997) . This 
observational breakthrough has, for the first time, provided 
both a precise localisation of a GRB position and broad- 
band spectral information. But, while the spectral fit by 
Wijers, Rees & Meszaros (1997) does agree quite well with 
the cosmological blastwave theory (Rees & Meszaros 1992), 
a coherent picture of the distance scale to GRBs has yet to 
emerge. Preliminary results as to the nature of the observed 
field surrounding the optical transient position have been 
inconclusive: a faint extended source near the optical tran- 
sient position may have faded from images taken at Keck at 
two different epochs (Metzger et al. 1997a; although see Fox 
et al. 1997). The natural assumption was that if GRBs are 
cosmological in origin then the optical transient would be 
associated with a galaxy; if the extension is indeed variable, 
then this is not the case. Furthermore, a spectrum of the 
optical transient and the extension taken with the Keck II 
telescope (Tonry et al. 1997) does not reveal any tell-tale 
features which would indicate the distance to GRB 970228. 



So far no extended emission has been detected from the re- 
cent burst GRB 970508 but Metzger et al. (1997b) have 
reported absorption lines in the spectra of the optical tran- 
sient putting the transient beyond z = 0.835. If confirmed, 
this would establish GRBs as cosmological. 

If GRBs are cosmological, what are the other likely 
observational consequences which would rule out the non- 
cosmological scenarios? One answer is that in most cosmo- 
logical models GRBs will have a parent galaxy which we 
can search for. Historically this has been difficult due to the 
scarcity of small error boxes: the Interplanetary Network 
(IPN) of satellites yielded only a few error boxes smaller 
then about 5 arcmin^ (Atteia et al. 1987). As the projected 
density of galaxies is reasonably high (eg. we expect about 
30 galaxies in 5 arcmin^ boxes at B < 24.0) there was no 
lack of potential hosts galaxies to the bursts; statistical stud- 
ies of GRB positions were thereby limited to searching for 
correlations on scale lengths greater than several arcmins 
(Kolatt & Piran 1996; Hurley et al. 1997). 

Kolatt & Piran 1996, for instance, do find evidence of 
a correlation of bursts with Abell clusters using the crude 
BATSE localisations (although see Hurley et al. 1997). Us- 
ing the better, but fewer, localisations from the IPN, there 
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have also been claims of a bright galaxy excess near bright 
bursts (Larson, McLean, & Becklin 1996) but this analysis 
uses Poissonian errors to quantify over- density, neglecting 
the fact that galaxies are clustered and counts far galaxies 
outside the strict limits of the IPN error boxes. Fenimore 
et al. (1993), examining 8 of the smallest IPN error boxes, 
conclude that host galaxies of GRBs must be fainter than an 
absolute visual magnitude M„ ~ —18. Recent Hubble Space 
Telescope (HST) observations of IPN boxes confirm this re- 
sult (Schaefer et al. 1997) and further constrain the host 
galaxies of GRBs to be fainter than absolute Bmag ^ —17.4. 
The conclusion that this requires the hosts of GRBs to be 
significantly fainter than L« is valid only if GRBs are stan- 
dard candles and if the faintest BATSE bursts have z <.l. 
Moreover, it does not allow for the possibility that a galaxy 
outside IPN error boxes could be the host to the GRBs. 

With the detection of the optical transient to GRB 
970228 we now have the beginning of a set of very high pre- 
cision positions and correlation analysis should thus begin to 
probe smaller scales at deeper magnitudes. Deep HST imag- 
ing reveals faint extended emission adjacent to the 970228 
optical transient point source. If further observations prove 
that this is a distant galaxy, then it must be a very good 
candidate for the host of the GRB. However, at a magni- 
tude of Mf606w — 25.5 (Tanvir & Johnson 1997) the spa- 
tial density of galaxies is high and the alignment could be 
a chance coincidence. Thus Galactic models are not ruled 
out and neither are cosmological models in which the pro- 
genitors of GRBs are ejected from their host or even only 
be loosely associated with the light density from galaxies. 
If GRBs arise from neutron star-neutron star (NS-NS) co- 
alescence, for instance, then a GRB may occur tens of kpc 
outside of its host galaxy (Sigurdsson, private communica- 
tion); at redshifts of z ~ 0.2, such positional displacements 
allow an offset of the optical transient from its host galaxy 
by as much as 20 arcsec. If GRBs do occur far outside their 
hosts then unambiguously identifying the parent galaxies 
will be difficult. However, to rule out the Galactic scenario, 
it would suffice to find a statistically significant over-density 
of potential host galaxies near GRB positions since bursts 
originating from within our Galaxy should be randomly dis- 
tributed with respect to underlying galaxy field. 

Clearly, to distinguish between the two distance scales, 
there must be a method which can both determine an over- 
density of potential hosts (if there is one) and accurately es- 
timate the significance of such an over-density. Ultimately, 
if the distance scale to GRBs is established to be cosmolog- 
ical, then the same analysis will help reveal how GRBs are 
distributed near galaxies. 

In section 2 we provide a framework for making quantifi- 
able statistical statements about the degree of over-density 
of hosts for a given number of GRB positions that accounts 
for galaxy clustering. In section 3 we describe the simulation 
of a galaxy distribution which is used in section 4 to pre- 
dict the number of GRB positions needed to find a signifi- 
cant correlation if GRBs are indeed associated with distant 
galaxies. In section 5 we show how the method developed 
herein may be extended to other types of data and hypothe- 
ses. In other words, should observations provide no evidence 
for a positive correlation, the test developed in section 2 can 
be used to accept the alternative hypothesis that GRBs are 
not correlated with potential hosts to faint limits (obviously. 



by rejection of the null hypothesis that there is a correlation 
under some given set of assumptions) . We consider this issue 
in Ught of results from HST images of GRB 970228. 



2 CONSTRUCTION OF A STATISTIC OF 
OVER-DENSITY 

2.1 What is a potential host? 

In an overwhelming number of viable cosmological models, a 
GRB originates from near or within a galaxy. In each model 
scenario, there is a three-dimensional probability, ^3(2?), of 
a GRB originating from a given galaxy at a co-moving po- 
sition, X, relative to its centre. In the simplest case, we take 
the probability of a GRB occurring to be constant within 
a sphere of radius I. The projected two-dimensional proba- 
bility distribution on the sky, for a galaxy at redshift z is 
thus; 



P2{i,e,z) 



where 



D 



IHo 



(1 



2c [(1 + 2) -(1 + ^)1/2] 



radians 



(1) 



(2) 



in the standard, qo = 1/2, A = 0, cosmology. Hereafter, 
the Hubble constant is taken as Hq — 70 km s~^ Mpc~^ 
(eg. Tanvir et al. 1995). For each burst, i, there is a to- 
tal probability normalisation, Ai{z), for each galaxy in the 
field. If the results of the log TV-log P studies are to be used 
as a guideline (as in Fenimore et al. 1993), then a bright 
GRB — such as the expected case for most BeppoSAX X-ray 
detections or small IPN error box localisations — is expected 
to come from a low-redshift galaxy and thus Ai{z) for that 
particular burst will peak at a low redshift. 

For each burst position i (which is observed in a given 
passband, and limiting magnitude) one constructs a prob- 
ability host map, PHMi(r), by adding the contribution of 
P2{r, 9, z)d9 from each galaxy where r is the maximum off- 
set scale of the map. A PHM gives the relative probability of 
a GRB occurring at any position on the sky assuming a par- 
ticular functional form of the offset probability, P2{r, 9, z). If 
GRBs are not correlated with galaxies then their positions 
should occur randomly on the PHM; conversely, if GRBs 
are correlated with galaxies then their positions will arise in 
regions of larger probability density on the PHM. 

Since we do not know a priori both the intrinsic offset 
scale I and the functional form of the offset probability dis- 
tribution, we must assume some map scale r and a form of 
the offset probability distribution [P2{r,6,z)]. In addition, 
construction of the PHM requires knowledge of the redshift 
of each galaxy; redshifts may be determined from photomet- 
ric colours, via their observed spectra, or simply estimated 
from their brightness using previous flux/redshift studies. 
Alternatively, a median redshift of galaxies in the field may 
be assigned to each galaxy. In the simulations in section 3, 
we show that, not surprisingly, a PHM is most sensitive to 
detecting correlations if r ~ We also find that assigning 
the median redshift to each galaxy in the map suffices in 
finding correlations (see fig. p|) at least for small intrinsic 
offset scales. Note that if the map scale r is sufficiently small 
(i.e. if GRBs are retained within a few kpc of a host galaxy). 
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the PHM reduces to a series of (5-functions. In section 6 we 
discuss what a reasonable value for I might be for different 
models. 

In what follows, let pij be the value of PHMi(r) for the 
jth randomly selected location on map i; by counting the 
frequency of pij for a large number of randomly selected 
points we construct the distribution fi{p), which is effec- 
tively the probability of finding p of a potential host on the 
sky at any point on the map. 



2.2 The Test for Over-densities of the True 
Sample 

For a given passband and limiting magnitude of an observa- 
tion of each position, there will be a corresponding number 
density, Ni, and two-point correlation function, LUi{0), that 
characterises the underlying distribution of galaxies from 
which fi{p) is found. In section 3 we simulate an underlying 
galaxy distribution that matches the observed properties. In 
practice, it may be better to estimate fi{p) directly from a 
number of control fields taken at similar airmass, galactic 
latitude, limiting magnitude, etc. 

The normalised cumulative probability function is given 



Hp) = 



(3) 



where pmin,i and Pmax,i are the extreme values of p on a 
given PHMi. If pi,obs is the observed value of p on the PHM 
of burst i then the set of numbers S — {-Fi(pi,obs)} will be 
evenly distributed in the interval [0,1] if the positions of the 
GRB ensemble are uncorrelated with that of the potential 
host galaxy distribution. 

The degree of deviation of S from an evenly distributed, 
cumulative set can be used to estimate the strength of 
the correlation between GRBs and observed potential host 
galaxies. Thus, this method allows one to combine an inho- 
mogeneous set of observations to test the hypothesis that 
GRBs are correlated with host galaxies. 

The significance of the discrepancy between the ob- 
served set S and the expected set is evaluated using the 
Wilcoxon statistic (see Noether 1967 and Efron & Petrosian 
1996). It is defined as: 



w = 



0.5 



(4) 



and compares the average value of S and the expected value 
of the cumulative probability distribution (= 0.5). The dis- 
tribution ofW values, ,f{W), for an uncorrelated set of GRB 
positions and potential hosts should be normally distributed 
with variance unity centred about zero. In section 3 we con- 
struct f{W) for a given limiting magnitude and number of 
optical transient source positions by a Monte Carlo simula- 
tion of many expected (uncorrelated) sets (see fig. [^). We 
test the null hypothesis that GRB optical transient source 
positions are uncorrelated with host galaxies by comparing 
the stochastic realization of Wohs with the expected distri- 
bution of f{W). Clearly if there is a positive correlation, 
then the M^-statistic will be negative (and the opposite is 
true for a anti-correlation). In this formalism, the null hy- 



pothesis can be rejected with a quantifiable confidence level 
by comparison of VKobs with f{W). 



3 PREDICTIONS 

In this section we construct simulations which will be used 
to assess the power of the above analysis under different 
assumptions about the true GRB distribution. Specifically 
we will use these simulations to ask how many observations 
would be required to show some correlation of GRBs with 
galaxies in various scenarios. 



3.1 Constructing a Probability Host Map (PHM) 

To estimate the expected significances, we construct a 
galaxy distribution which artificially replicates the observed 
galaxy number density, two-point correlation function, and 
redshift distribution for a given limiting magnitude. Several 
galaxy distribution studies have provided the necessary data 
to construct such a distribution but here we limit ourselves 
to the correlation/redshift study of Roche et al. (1996). For 
a given limiting magnitude m in the R- and B-band, they 
provide a redshift distribution of galaxies {fm{z)), number 
density (Nm), and fit to a two-point correlation function 
{ujrn{0)). The observations were made at high galactic lat- 
itude (6 ~ 85°) and with little extinction {Eb — 0.15 and 
Er ~ 0.08). This simulation is obviously extendible using 
other galaxy distribution studies which provide similar data 
in other passbands and/or limiting magnitudes. 



3.2 Simulating a PHM from galaxy spatial and 
redshift distribution 

As mentioned earlier, it may be possible to construct a PHM 
directly from control fields taken under similar conditions 
as the optical transient field. Here, however, we generate a 
PHM artificially. To generate an artificial galaxy distribu- 
tion, we employ the truncated hierarchy method (Soneria & 
Peebles 1978) which has been adapted to 2-D projections 
(Infante 1994). For a given burst position i we vary the in- 
trinsic scaling of the algorithm (see Infante 1994) to fit the 
expected spatial distribution of underlying galaxies (quan- 
tified by the spatial density distribution, Ni and the ampli- 
tude and slope, a, of the two-point correlation function). For 
each galaxy in the distribution, we construct a probability 
profile from eq. |l] replacing I with r and thereby construct 
a PHM(r). 

The construction of the galaxy distribution and PHM 
requires a few assumptions. 

(i) We take the slope of the two-point correlation func- 
tion to be a = —0.8 which, based on numerous galaxy dis- 
tribution studies (eg. Roche et al. 1996), seems reasonable. 
In our simulation using the truncated hierarchy method, for 
B < 24.0, we are able to reproduce the two-point correlation 
function and number density {N — 20730 gals/deg^) that 
match that found in Roche et al. 1996. Bright GRBs are ex- 
pected to come from moderate redshifts. We have therefore 
chosen a limiting magnitude of B < 24.0 which corresponds 
to a median redshift of 2 ~ 0.6. Clearly the analysis may be 
extended to fainter magnitudes in which case the expected 
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Figure 1. Simulated probability host map (PHM) from galaxy fields {B < 24.0) for a) r = 50 kpc, r = 200 kpc, and, c) r = 500 kpc. 
Regions of higher potential host density are darker and the size of the circular region around each galaxy is related to its redshift by 
eq. [0. Each field, comprising about 2300 galaxies, is a 20 arcmin X 20 arcmin portion of a large PHM. The simulated clustering of 
galaxies evident gives a larger spread in /(p) than that derived from a randomly distributed set of galaxy positions. 



number of burst positions required to find a correlation will 
increase. 

(ii) We assign a redshift to each galaxy in the generated 
distribution. Roche et al. (1996) provide an estimated galaxy 
redshift distribution as a function of limiting magnitude. 
For each cluster in an artificially generated spatial distri- 
bution, we then assign a similar redshift (±5 percent) to its 
constituent galaxies based on a probability distribution con- 
structed from the Roche et al. (1996) data. This accounts for 
both for the observed distribution of redshifts and the fact 
that galaxies within the same cluster have similar redshifts. 

(iii) We take Ai{z) equal to unity. Note that any non- 
uniform dependency of Ai{z) upon redshift will reduce the 
spatial density of potential hosts and thus increase the de- 
tectability of a correlation. Although the distances to GRB 
hosts should be in accord with distance implied by the log N- 
log P, it is not clear that the log N-log P gives a meaningful 
constraint on the functional form of Ai{z) since a standard 
luminosity of GRBs appears to be unconstrained (Loredo & 
Wasserman 1996; Horvath, Meszaros, & Meszaros 1996). 

A 20' X 20' portion of three PHMs is depicted in fig- 
ure (jl|) for various map scale lengths, r. The clustering of 
galaxies is evident. 



4 SIMULATION RESULTS 

At a given limiting magnitude, bandpass, and number of op- 
tical transient positions, A^'burst, we find the Wilcoxon distri- 
bution fiW) as prescribed in section II. We then simulate 
an over-density of host galaxies near GRBs by choosing a 
random galaxy on the sky and placing the optical transient 
position at an offset prescribed by eq. [Q]. The stochastic 
realization of H^obs vs. /(W) is used to determine the signif- 
icance of a correlation. Figure (^ shows an example com- 
parison of Wobs (dashed vertical line) vs. /(W). For a given 
limiting magnitude, map scale r, and intrinsic scale I, we 
predict the probability that the stochastic realization of W 
is consistent with the null hypothesis, 

Hq: GRB positions are uncorrelated with potential hosts. 

For a magnitude B < 24.0 and r = 30 kpc, figure (S) shows 
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Figure 2. The expected distribution of f{W) for an uncorrelated 
set of 6 GRB positions for I = r = 100 kpc. Also, shown with a 
vertical dashed line, is the predicted W determined from a corre- 
lated set of 6 GRB positions. From the distribution, we compute 
that such a value {W ~ —2) or smaller could occur by random 
chance in only 3.1 percent of cases if the GRB positions are not 
correlated with potential hosts. Thus, the null hypothesis may be 
rejected at greater than 95 percent confidence. 

the confidence in rejecting Ho as a function of the number 
of GRB positions and I. Notice that if I is less than ~ 50 
kpc, then only a few GRB positions are required to rule out 
Ho and conclude that GRBs are correlated with host galax- 
ies. This is not a surprising result given that the typical 
projected distance at z ~ 0.3 for I = 30 kpc is only a few 
arcsec: the expected value of PHM(r — 30) is close to zero 
and thus any over-density is rare. The correlation analysis 
is sensitive to the chosen map scale r and whether knowl- 
edge of galaxy redshifts are included. Figure (^ shows the 
expected number of positions needed to find a 95 percent 
confident over-density (if there is one) as a function of I and 
r where the PHMs were constructed using only galaxy po- 
sitions and the same (median) redshift {z ~ 0.61) for each 
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Figure 3. Null hypothesis rejection confidence as a function of the number of GRB positions. The null hypothesis tested is that GRBs 
are uncorrelated with the potential host distribution. The PHM is made with r = 30 kpc and taking either a) every galaxy is located 
at the median redshift (2 ~ 0.61) of B < 24.0 galaxies (Roche et al. 1996) or b) the simulated redshift for each galaxy. The predicted 
curves for I = 30 kpc (solid), I = 50 kpc (dot), I = 100 kpc (dash), I = 200 kpc (dot-dash), and I = 500 kpc (dots-dash) are depicted; the 
small deviations from a smooth curve are the result of a finite number of trials in the Monte Carlo simulations. For I < 50 kpc, only a 
few burst positions are required to find a 99 percent confident correlation of GRBs with potential hosts. Note that with increasing values 
of / it becomes more difficult to detect over-densities. However, other chosen PHM scales r are more sensitive to correlations at larger I 
(see fig. [^). The predictions do not greatly differ for each type of map. Thus, precise knowledge of the redshifts of galaxies will not be 
necessary to construct a PHM. 



galaxy. Small values of r produce PHMs that are most sen- 
sitive to a correlation if I is also small; likewise, correlations 
are most readily found if / is large by a choice of a large value 
of r in making the map. As evidenced by the small number 
of bursts required to make a significant correlation, the red- 
shift of each galaxy, inclusion of which strengthens the test 
for correlations, need not be precisely known in making the 
PHM. The simulation predictions presented in fig. are 
made from PHMs where the same single redshift is assigned 
to all galaxies. Figure ^ illustrates the minor strengthening 
of the correlation test where knowledge of galaxy redshifts 
is included in construction of the PHM. 



20 r 
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5 EXTENSION TO OTHER DATA SETS AND 
HYPOTHESES 

5.1 Using IPN localisations 

A significant correlation of GRBs with observed galaxies 
may be found without further precise locations that opti- 
cal transients would provide. Instead, the method described 
above may be used on observed fields near existing IPN lo- 
calisations. In such a case, fi{p) is estimated by randomly 
selecting positions on a PHM and taking p to be the aver- 
age density within an error box of identical shape and size 
as IPN localisation i. The observed value of pi is then com- 
pared against fi{p) as above. 



50 100 150 200 

k'ap Scale r (kpc) 

Figure 4. Number of GRB positions needed for 95 percent con- 
fident correlation detection as a function of map scale r. The 
intrinsic scale of GRBs in the simulation are I = 30 kpc (solid 
line), I = 50 kpc (dotted), I = 100 kpc (dashed), and I = 200 
kpc (dot— dashed). Notice that a correlation is most readily found 
with r ^ I. 



5.2 Are GRBs Uncorrelated with Observable 
Galaxies? 

If the extension near the position of GRB 970228 is not a 
galaxy there are no hosts at the faint EST limit (AdprnKW < 
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Table 1. Gamma-ray burst 970228 constraints on the minimum offset of GRBs from their host galaxy. The data found is from cleaned 
HST PCI images taken on March 26 and April 7 (Tanvir & Johnson 1997). Instrumental magnitudes are derived in two passbands and 
errors are found from systematic (~ 0.1) and photometric uncertainties. Zero-point calibration is taken from Holtzman et al. (1995). 
Redshifts are estimated from an F606W magnitude comparison with galaxies in the Hubble deep field (Cowie 1997). The redshift at 
which 90 percent of the Cowie galaxies of similar magnitudes have larger photometric redshifts is denoted by 290 (similarly, 2:50 is the 
median redshift of galaxies with similar magnitudes). This corresponds to a constraint on the minimum offset Iqq (assuming go = 0.5) of 
GRBs from their host galaxies. 



Galaxy 
Number 


Position (J2000) 
RA (hms) DEC (deg) 


D 
(arsec) 


Ma; 

F606W 


gnitude 

F606W - F814W 


Redshift 
Z50 (Z90) 


'50 ('90) 
{h~Q kpc) 


Galaxy Property 


1 


5:01:46.50 


11:46:50.2 


4.1 


26.3 ±0.4 


0.9 


1.34 (0.25) 


25 (14) 




2 


5:01:46.20 


11:46:45.4 


10.8 


21.8 ±0.1 


1.1 


0.31 (0.09) 


43 (17) 


spiral, star-forming? 


3 


5:01:46.56 


11:46:46.1 


7.7 


23.7 ±0.1 


0.8 


0.70 (0.20) 


44 (23) 




4 


5:01:47.17 


11:46:49.3 


8.8 


25.5 ±0.21 


0.8 


0.88 (0.62) 


53 (48) 




5 


5:01:47.15 


11:46:54.4 


5.2 


25.2 ±0.23 


1.0 


0.87 (0.25) 


31 (18) 


emission-line galaxy 














z = 0.638" 


28.9 




6 


5:01:46.94 


11:47:04.1 


11.4 


25.4 ±0.30 


0.5 


0.84 (0.61) 


67 (62) 





" Spectroscopic redshift from Tonry et al. 1997. 



26.5) for a GRB/host-galaxy offset of I < Uh^Q kpc (see ta- 
ble 1). If we take the less faint limits of B < 24.0 commen- 
surate with simulations, then the somewhat stronger con- 
straint of Zmin > i7h^Q kpc comes from galaxy 2 in table 
1. As the redshifts of the galaxies in the HST fields are yet 
to be determined spectroscopically, these offset lower limits 
should be considered approximate. 

While this one optical transient position does not pro- 
vide a strong constraint, it does seem that either it has a 
surprisingly faint host or GRBs occur at large offsets from 
their host galaxy. If the extension is indeed a galaxy, then it 
is so faint that there is a non- negligible chance (few percent) 
that the optical transient position is coincident with it by 
random chance. 

The absence of a detectable galaxy down to faint mag- 
nitudes {Mf606w < 25) near the optical transient position 
of GRB 970228 gives sufficient impetus to the formulation of 
the alternative null-hypothesis that GRB positions are cor- 
related with potential hosts. A significant rejection of this 
hypothesis will place important constraints on most cosmo- 
logical models. Within such a hypothesis, f{p) is constructed 
by choosing positions near galaxies (rather than at random) 
from a PHM. By examination of the HST fields taken on 
26 March and 7 April, already one can rule out Ho for 
I < 17h:^o kpc for B < 24 (galaxy 2; table 1). Obviously, 
GRBs could originate from low-surface brightness galaxies 
or at large redshifts, in which case the detectability of a 
correlation would be difiicult. 



6 DISCUSSION 

Regardless of the interpretation of the emission near the op- 
tical transient position, it is clear that GRB 970228 did not 
originate from a bright, nearby galaxy (see table 1) con- 
firming the systematic findings of Schaefer et al. (1997) and 
Fenimore et al. (1993). It should be emphasised that the "no- 
host" conclusion of Schaefer et al. (1997) and Fenimore et 
al. (1993) is based strongly on the assumption that the GRB 
population follows a non-evolving, standard candle luminos- 
ity in as much as it is argued that bright bursts should come 



from nearby galaxies. But from an observational standpoint, 
Loredo & Wasserman (1996) have concluded that the width 
of the luminosity function of GRBs is poorly constrained 
and Horvath, Meszaros, & Meszaros (1996) find that a large 
range of luminosities acceptably fit the log A^-log P distri- 
bution. Moreover, theories of cosmological shocks and blast- 
waves predict a GRB luminosity which can vary greatly 
from burst to burst (eg. Fenimore, Madras, & Nayakskin 
1996); instead total energy is more likely to be standard (see 
Bloom, Fenimore, & in 't Zand 1997). Thus, robust conclu- 
sions about (non-) correlations of hosts with GRBs should 
not rely on the above assumptions. 

The specific correlation hypothesis that is tested is 
wholly determined by the functional form and dependences 
in eq. . The basic hypothesis we have tested in these sim- 
ulations is that "the progenitors of GRBs originate (and are 
possibly ejected) from observable galaxies." If mass in the 
Universe scales with light, this can also be seen as a test 
that the rate density of GRBs is highest at the peaks of 
the matter distribution. Of course while this hypothesis is 
somewhat model dependent, it effectively introduces only 
one free parameter, namely I (and to a lesser extent P2). 

We conservatively estimate the maximum intrinsic off- 
set scale to be of order I = 400 kpc. This corresponds to a 
constant progenitor outward velocity «pro = 400 km/s and 
a collision time scale of r = 10® years. If all GRBs occur 
at such large distances from their host galaxy, it would be 
difficult to find a significant correlation with a small collec- 
tion of GRB positions. However, simulations have shown for 
NS-NS models that / may be an order of magnitude smaller 
and P^il) is likely to be very peaked near the centre for most 
galaxies, so that the effective I might be significantly lower 
than 400 kpc (Lipunov et al. 1995). Therefore, whatever the 
true / is, we believe that the assumed P2{r — I), is quite 
conservative by over-predicting the true extension of GRBs 
from their hosts. 

Since / is not known it is impossible to know a priori 
which map scale r will optimise the detection of a correla- 
tion. The predictions of a 95 percent confident correlation af- 
ter about six optical transient detections is based on a fixed 
r = 60 kpc and assuming I < 100 kpc (see fig. [4]). Although 
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a detected correlation thusly will be adequate for resolving 
the distance scale to GRBs debate, ideally one would like to 
determine I in addition. As a correlation search over a range 
of map scales r is likely to increase the random chance of 
finding a correlation, it will be necessary to determine I a 
posteriori via Bayesian analysis. 

One limitation of the statistical method presented is 
that the angular size of the PHM must be significantly larger 
than the projected angular offset D of GRBs from their host. 
As most cosmological models would have GRBs originate 
at distances significantly less than 1 Mpc from their host 
galaxy, this corresponds to a minimum PHM size of about 
30 arcmin x 30 arcmin. The correlation analysis hinges on 
an accurate modelling of the expected probability density 
distribution (/(p)) and thus control fields must be selected 
quite carefully in order not to bias the PHM. We will exam- 
ine the salient issues of estimating /(p) in a future paper, 
which presents the results of a deep optical IPN error box 
correlation analysis. 



7 CONCLUSIONS 

We have presented a formalism which may be used to quan- 
tify the degree of over-density of potential host galaxies near 
GRB positions; if GRBs are found to be cosmological in ori- 
gin this formalism may also be used to determine the intrin- 
sic offset of GRBs from their hosts. The broad hypothesis 
we have chosen to address is that GRBs originate within / 
kpc of their host galaxy. In our simulations, we find that for 
I <^ 100 kpc, encompassing most predictions of current cos- 
mological models, a 95 percent significant correlation may be 
found with about six positions. This prediction corresponds 
to simulations of galaxy fields of B < 24.0 where the median 
redshift of galaxies is z ~ 0.6. If GRBs are correlated with 
galaxies at comparable redshifts, such correlation may found 
within a year if BeppoSAX detections lead us to a handful 
of optical transient localisations. 

However, the extension near the optical transient of 
GRB 970228 is very faint {Mpeoew ^ 25.5 implying a red- 
shift 2 ~ 1), which suggests that GRBs may only correlate 
with very distant galaxies. In this case, correlations will be 
detectable only if I is sufficiently small. If the fading ex- 
tension near the optical transient of GRB 970228 is not a 
galaxy, we find that the progenitor system must have trav- 
elled a minimum distance Imin ^ i4:h^Q kpc from its host 
galaxy before producing a GRB. Assuming that galaxy de- 
tection in the HST exposures is complete to distances of 
the bright burst, this effectively requires GRBs to originate 
well outside the nucleus of their host galaxies. In any event, 
the absence of a bright host galaxy indicates that the non- 
evolving standard luminosity derived from the log A'^-log P 
distribution inadequately describes the true distribution of 
cosmological GRBs. 
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